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We show that the gamma ray spectrum observed with the HESS array of Cherenkov telescopes 
coming from the Galactic Center (GC) region and identified with the source HESS J1745-290, is 
well fitted by the secondary photons coming from dark matter (DM) annihilation over a diffuse 
power-law background. The amount of photons and morphology of the signal localized within a 
region of few parsecs, require compressed DM profiles as those resulting from baryonic contraction, 
which offer ~ 10 3 enhancements in the signal over DM alone simulations. The fitted background 
from HESS data is consistent with recent Fermi-LAT observations of the same region. 



Observations of very high energy (VHE) 7-rays coming 
from the Galactic Center (GC) have been reported by dif- 
ferent collaborations such as CANGAROO Q], VERITAS 
0, HESS I3IH, MAGIC [5] and Fermi-LAT [S[7j. In this 
work, we will focus on the data collected by the HESS 
collaboration during the years 2004, 2005, and 2006 asso- 
ciated with the HESS J1745-290 source [J. The absence 
of variability in the TeV data suggests that the emis- 
sion mechanism and emission regions differ from those 
invoked in the variable IR and X-ray emission [5]. Im- 
portant deviations from a power law spectrum has been 
already proved, and a cut-off at several tens of TeVs is a 
remarkable feature in the data. The angular distribution 
of the VHE 7-ray emission of HESS J1745-290 shows the 
presence of an adjunctive diffuse 7-ray emission compo- 
nent, but the significance of the signal reduces to few 
tenths of degree in any case [4]. 

The fundamental nature of this source is still unclear. 
These gamma rays could have been originated by parti- 
cle propagation [7J |9] in the neighborhood of the Sgr A 
East supernova remnant and the supermassive black hole 
Sgr A, both located at the central region of our galaxy 
[TUl [TT] . In this letter, we will analyze the possibility 
of explaining the spectral features of the signal with the 
photons produced by the annihilation of DM particles. 
This interpretation has been widely discussed in the lit- 
erature from the very early days of the publication of 
the observed data by the above collaborations [12j [13] . 
It was concluded that the spectral features of the data 
disfavored the DM origin |13j . However, in our study, we 
will assume that the DM signal is not the only contribu- 
tion, but it is complemented by a background which is 
well motivated by radiative processes generated by par- 
ticle acceleration in the vicinity of Sgr A East supernova 
and the supermassive black hole: 
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In order to simplify the parameters in our fits, we will 
assume just a simple power-law for the contribution not 
associated with DM: 
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This assumption is also experimentally motivated by 
Fermi-LAT data corresponding to 25 months of observa- 
tions of the source IFGL J1745. 6-2900, that is claimed to 
coincide spatially with the HESS J1745-290 source [Ti] . 
These data have been shown to be well fitted by a bro- 
ken power-law, that for E > 2 GeV is consistent with the 
spectral index T = 2.68 ± 0.05 {x 2 /dof = 0.81) [7j. In 
any case, we will not assume any prior on the parame- 
ters for the background, and either B and T will be fitted 
from HESS data in our analysis. 

The differential gamma-ray flux coming from DM par- 
ticles can be written in general as: 
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where a = 2 takes into account the gamma rays coming 
from DM annihilation (we are assuming that the DM 
particle is its own antiparticle) , whereas a = 1 accounts 
for the photons generated by the possible DM decays. To 
determine the Galactic signal, it is necessary to compute 
the astrophysical factors (J( a )) in the direction given by 
the '5 angle, defined as the one between the direction of 
the Galactic center and the line of observation, given by: 
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where I is the distance from the Sun to any point in the 
halo. The radial distance r is measured from the GC, and 
is related to I by r 2 — I 2 + Dq — 2D Q l cos where D Q ~ 
8.5 kpc is the distance from the Sun to the center of the 
Galaxy. The distance from the Sun to the edge of the halo 

in the direction 9 is l max — D cos 9 + ^Jr 2 — Dq sin9. 

The astrophysical factor is proportional to p 2 when ac- 
counting for DM annihilation channels (whereas it is just 
proportional to p when computing photons from DM de- 
cays). 

The photon flux is maximized in the direction of the 
GC, and must be averaged over the solid angle of the 
detector. For detectors with sensitivities in the TeV 
regime, the solid angles are typically of order Afi = 
2tt(1 - cos*) ~ 10~ 5 , as it is the case for the HESS 
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FIG. 1: Best fit to the HESS J1745-290 collection of data (years 
2004,2005, and 2006 |4|) in the case that the DM contribution came 
entirely from annihilation into uu quarks. The full line shows the 
total fitting function (Eq. ([!})■ The dotted line is the fitted power- 
law background given by Eq. The dot-dashed line corresponds 
to the DM annihilation contribution with resolution uncertainty 
(R.U.) of 15%, typical from HESS (the DM contribution without 
R.U. is shown by the dotted line for reference). It is remarkable to 
note that the cut-off in the spectrum characteristic of this annihi- 
lation channel coincides with data. The parameters of the fit are 
reported on Table [I] 



Cherenkov telescopes array. The dark halo is usually 
modeled by the NFW profile [T5], that is in good agree- 
ment with non-baryonic cold DM simulations. However, 
it has been claimed [TH [17] that when baryonic gas is 
taken into account, it falls to the central region, modi- 
fying the gravitational potential and increasing the DM 
density in the center (see however [TS]). This fact has 
two important consequences for our analysis. On the 
one hand, the central region accessible to gamma ray de- 
tection is compressed to few tenths of a degree; on the 
other hand, the DM annihilating fluxes are enhanced by 
a factor ~ 10 3 over the classical NFW profile 17 J. The 
interpretation of the HESS data from DM annihilation is 
in good agreement with these types of compressed dark 
halos. 

On the other hand, the rest of the computation de- 
pends on the nature of the DM particle. If it is meta- 
stable, photons can be produced by its decay. In this case 
(a = 1): C| = Tj is the decay width into SM particles 
(labeled by the subindex i). This possibility is not com- 
mon to all DM candidates. Much more general is the DM 
annihilation since DM particles are their own antiparticle 
in the most part of DM models. The case a — 2 takes into 
account the gamma rays coming from DM annihilation 
since Q = {(Jiv) are the thermal averaged annihilation 
cross-sections of two DM particles into SM particles (also 
labeled by the subindex i). M is the mass of the DM par- 
ticle, and the number of photons produced in each anni- 
hilating or decaying channel dN^ / dE , involves decays 
and/or hadronization of unstable products such as quarks 
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FIG. 2: Same as in Fig. [l]but with annihilating DM into 
gauge bosons. Electroweak annihilation channels are softer than 
hadronic ones. The large uncertainties for data over ~ 10 TeV do 
not allow to discriminate among these two types of spectra. Both 
of them are in good agreement with HESS observations. 
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FIG. 3: Same as in Fig. [l]but with annihilating DM into t+t~. 
The poor quality of the fit is evident and common for all the lep- 
tonic channels. 



and leptons. Because of the non-perturbative QCD ef- 
fects, the calculation of dN^ jdE requires Monte Carlo 
events generators such as PYTHIA [Tj5]- However, the 
fact that simulations have to be performed for fixed DM 
mass implies that we cannot obtain explicit M depen- 
dence for the photon spectra. In order to overcome this 
limitation, different fitting functions have been obtained 
in |20j for a wide range of masses which will allow us to 
include the mass M as an additional parameter in the 
fits to HESS data. Those fitting functions are shown to 
depend on several parameters whose scaling behaviour 
with M were also obtained in [50]. This fact has allowed 
us to extrapolate the simulated spectra to the very high 
energies required in this analysis. 

In this letter, we will focus on gamma rays coming 
from external bremsstrahlung and fragmentation of SM 
particle-antiparticle pairs produced by DM annihilation. 
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TABLE I: In this table, the four parameters of the annihilating DM into a single channel fit are presented: The fitted value of 
the mass (TeV) of the annihilating WIMPs, normalization factor of the signal A (10 -7 cm"^ -1 ' 2 ), normalitazion factor of the 
gamma ray diffuse emission background B (10 -4 GeV _1/ ' 2 cm _1 s _1/ ' 2 ), and spectral index V of the same background. The \ 2 
per degree of freedom (dof), and the value of its variation A\ 2 = X 2 ~ Xdd with respect to the best one (Xdd) i s a l so provided 
(In the case of four parameters, values of A\ 2 = 4.72,9.70 and 13.3 correspond to 68.3% , 95.4% and 99.0% confidence level, 
respectively). Finally, the astrophysical factor is computed by assuming (av) = 3 • 10 -26 cm 3 s _1 , and presented in units of the 
the astrophysical factor associated with a NFW profile: 6 = (J(2))/<^) ?W ), where (J^ w ) ~ 280 • 10 23 GeV 2 cm -5 . 



We will ignore DM decays, the possible production of 
monoenergetic photons, n-body annhilitations (with n > 
2), or photons produced from internal bremsstrahlung, 
that are model dependent. In particular and in order to 
simplify the discussion and provide useful information for 
a general analysis, we will consider DM annihilation into 
each single channel of SM particle-antiparticle pairs, i.e. 
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is a new constant that will be fitted together with the DM 
particle mass M , and the background parameters B and 
r. We assume a typical experimental resolution of 15% 
(AE/E ~ 0.15) and a perfect detector efficiency. The re- 
sults of these 4-parameters fits are summarized on Table 
|l) We can see that the best fit is provided by the dd chan- 
nel with x 2 /dof = 0.74 for a total of 24 dof. In any case, 
other hadronic channels such as uu (see Fig. [I]) or ss, 
also provide very good fits within la. In the same way, 
softer spectra as the one provided by ZZ, W + W~ (see 
Fig. [5]) or tt channels are consistent with data without 
statistical significance difference. On the contrary, lep- 
tonic channels (not only e + e~, or but also t + t~, 

Fig. pjp, cc and bb channels are ruled out with more than 
99% confidence level when compared to the best channel. 

It is interesting to note that taking into account all 
the channels that provide a good fit, the DM mass is 
constrained to 15 TeV < M < 110 TeV within 2a. The 
lighter values are consistent with hadronic annihilations 
(uu) and the heavier ones with the annihilation in tt, that 
is more similar to electroweak channels. On the other 
hand, at the same 95% confidence level, the allowed range 



for the spectral index of the diffuse background is 2.4 < 
r < 3.7. In this case, the lower values are consistent with 
all the allowed channels, but the higher values are only 
accessible to the light quark channels. In any case, as 
shown in Table I, all the channels that provide a good fit 
to HESS data are also consistent with the spectral index 
observed by Fermi-LAT data: T — 2.68±0.05 for energies 
between 2 GeV < E < 100 GeV [7] as we have already 
discussed. 

If we know the value of the annihilation cross-section, 
it is possible to obtain the astrophysical factor. In 
Table |TJ (J(2)) is given for each channel by assuming 
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3-10 cixrs and presented in units of 



the astrophysical factor associated with a NFW profile: 
b = (J(2))/(J(^ W )- We see that these results are consis- 
tent with b ~ 10 3 , that are predicted by simulating the 
baryonic dissipation effect in the DM halo [TT] as we have 
commented above. 

For a particular DM candidate, several channels can 
actually contribute. In this situation, the above analy- 
sis can be helpful as a guide but not determinant, since 
a combination of channels can work in a different way. 
For example, an interesting DM candidate which could 
have high enough mass and account for the right amount 
of DM in the form of a thermal relic, is the branon 
|21j . associated with brane fluctuations in brane- world 
models. For masses over 1 TeV, the main contribution 
to the photon spectra comes from branons annihilating 
into gauge bosons ZZ, and W + W~ . They produce ap- 
proximately the same amount of Z , W + and W~ since 
(a w+w -v) ~ 2(a zz v) ~ A/ 6 /(8tt 2 / 8 ), where / is the 
brane tension scale (and we are assuming only one bra- 
non specie 21). As the rest of channels can be neglected, 
and the W and Z channels produce very similar pho- 
ton fluxes, the results for the fit with the branon model 
lead to very similar results to those obtained by consid- 
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TABLE II: Best fit parameters for branon annihilation. The dominant channels contributing are W + W and ZZ. In the 
second column C 2 = Afi ( J( a )}An/(8-7rM 2 ). b is computed with the model fitted cross-section to WMAP data [22]: (av) = 
(1.14 ±0.19) ■ 10" 26 cm 3 s"\ Rest of parameters are as in Table I. 



ering both channels individually (see Table II). In this 
case, as there is a particle model behind, we can deduce 
the coupling that leads to the DM abundance consis- 
tent with WMAP observations [22]: / = 27.5 ± 2.4 TeV, 
and we can compute the thermal averaged cross-section: 
(av) = Y,i=w,z( (T i v ) = (1-14±0.19)-10" 26 cmV 1 , that 
agrees with the expected order of magnitude. 

In this work, we have analyzed the possibility of ex- 
plaining the gamma ray data observed by HESS from the 
central part of the galaxy by being partially produced 
by DM annihilation. We have proved that even single 
channel annihilations provide good fits if the DM signal 
is complemented with a diffuse background compatible 
with Fermi LAT observations. The morphology of the 
signal is consistent with dark halos compressed by tak- 
ing into account baryonic dissipation [T7] . The DM 
particle that may have originated these data needs to be 



heavier than ~ 10 TeV. This makes extremely difficult 
that these particles could be observed in direct detection 
experiments or produced in particle colliders [22] • In this 
sense, the analysis of other cosmic rays [24] (not only pho- 
tons, but also antiprotons, positrons, neutrinos,...) from 
the GC and from other astrophysical objects is funda- 
mental to cross check the hypotheses considered in this 
work. 
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